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SOLAR RADIATION AT THE EARTH'S SURFACE: ITS CALCULATION 

AND INFERENCE FROM SATELLITE IMAGERY 

Fred L. Bartman 
University of Michigan 
Ann Arbor, Michigan 


INTRODUCTION: PHYSICAL AND EMPIRICAL MODELS 

Solar radiation incident at the surface has not been measured at all loca- 
tions on the Earth’s surface. A valuable supplement to measured values are 
calculated values. This paper discusses a few methods for calculation of the 
insolation on a horizontal surface and refers to some of the many papers in the 
literature dealing with this topic. 

A valuable reference for this paper has been the notes for the intensive 
short course entitled "Solar Energy Measurements and Instrumentation" which was 
presented as part of the University of Michigan Engineering Summer Conference 
Program on July 9-10, 1979. The notes were written by staff members in the 
Department of Atmospheric and Oceanic Science. Of particularly great help in 
the preparation of this paper is Chapter 8 of those notes, "Modelling of 
Received Solar Radiation," by W. R. Kuhn. 

The various calculations can be broadly grouped into two categories, physi- 
cal models and empirical models. Physical models use the basic physical princi- 
ples of radiative transfer. They require the input of those quantities involved 
in the attenuation of the solar beam as it passes through the atmosphere. Data 
on these factors, i.e., the vertical distribution of water vapor and ozone, the 
aerosol vertical distribution, particle size distribution, and index of refrac- 
tion, the physical characteristics of clouds, and the reflectivity of the 
Earth's surface, are not always available and are highly variable. In addition, 
the calculations are quite elaborate. Thus the other type of model, the empiri- 
cal one, is often desirable. Empirical models are generally developed by using 
a known data set to calculate a set of regression coefficients for some set of 
parameters, such as the solar zenith angle, the duration of sunshine, or the 
amount of opaque cloud cover. Then the model can be used ,to predict the inci- 
dent solar radiation for other values of the parameters. These models are easy 
to use but must be applied with care. 

The advent of satellite remote sensing of clouds and the surface make the 
use of empirical models employing satellite images possible. Such a technique 
is also discussed. 

In all the following sections, the surface receiving solar radiation is 
assumed to be horizontal. The case of solar radiation incident on an inclined 
plane is not discussed here. 
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PHYSICAL MODELS 


Under clear sky conditions, we can write an equation for the flux (F) of 
radiation arriving at the surface as 


F = I cos e + D (1) 

2 

where F is the global solar radiation (W/m ) , I is the direct beam of solar 
radiation, D is the diffuse component, and 0 is the solar zenith angle. 


Spectral Calculations 

In terms of spectral components, I and D above can be written as 


I = I exp(-Tj^ sec 0) dA 

Jo 

rco 

D = I dA 

Jo 


( 2 a) 


( 2 b) 


where spectral solar irradiance at the top of the atmosphere, 

is the spectral extinction optical depth, and A is wavelength (ym) . 

The extinction is due primarily to absorption by ozone (O 3 ) and water 
vapor (H 2 O) , to scattering by air molecules, i.e., Rayleigh scattering (R) , and 
to absorption and scattering by aerosols (A). Thus, 




(3) 


Leckner (1978) presented expressions for these extinction coefficients. The 
Rayleigh scattering extinction coefficient is given by 


T^^(R) = 0.008735A (4) 

Po 

where p is the surface pressure (mbar) and p^ is 1013 mbar. For ozone, 

X;^{ 03 ) = -k(A) SL (5) 

where k(A) are spectral absorption coefficients (cm ^) given in table I and 
£ is the total ozone amount in atm cm. 
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The optical depth for water vapor can be written as (see also McClatchey 
et al. , 1971) 


T^(H20) 


0.3 k^(A) 


[l + 25.35 k„(A) X^]°' 


45 


( 6 ) 


where effective absorption coefficients of water vapor (cm^-g~^) , 

is the total precipitable water (g/cm^) , and m = sec 0 , the optical air 
mass. The values of k^(A) are given in table II. 


The aerosol optical depth can be quite accurately written as 


Tj(S) 



Q, (z) Tr[r (z) ] dz 


(7) 


where n(z) is the aerosol number density at height z, mean 

spectral extinction efficiency for the aerosols at altitude z, and r(z) is 
the mean aerosol radius at altitude z. Since the vertical variations of aero- 
sol concentration, composition, size distribution, and index of refraction are 
difficult to know, it is best to use the simple power law expression for aerosol 
attenuation 


T^(A) = -6A 


( 8 ) 


where 0 is the turbidity coefficient and a is the wavelength coefficient. 

Measurements of a indicate a wide range of possible values. However, 
most values occur in the range from 0.8 to 2, with an average value of 1.3 being 
considered most reasonable. The turbidity coefficient can be measured with 
reasonable accuracy. Typical values are shown in table III. 

The spectral diffuse component of solar radiation can be written as 

D;^ = 0.5^I;^(0) exp|-|x;^ (O 3 ) + (H20)j m| - I;^^ COS 6 (9) 

Leckner (1978) has made calculations of this kind and compared them with experi- 
mental measurements in a case where 3 was known to lie in the range from 0.05 
to 0.1. The results are shown in figures 1 and 2. 


Calculations for the Total Integrated Solar Radiation 

The previous section indicates that the spectral calculation for solar 
radiation at the surface is fairly, complex. Some simplification is achieved 
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when the integrated values are calculated with wavelength- averaged values of 
the parameters involved and with aerosol effects neglected. 

An effective mean value of Rayleigh optical depth "TqCR) is given in 
table IV as a function of zenith angle. The tabulated values refer to a verti- 
cal path only. Thus the direct solar radiation reaching the surface if the only 
extinction process were Rayleigh scattering would be 


F(R) 


cos 0 exp 



(R) 


Tan 


(10) 


where Iq is the solar constant, 6 is the zenith angle, and p is the sur- 
face pressure in mbar. The values of table IV can be calculated from 


Te(R) 


1 

sec 0 



^ OO 

I lx<0) dX 

Jo 


l^iO) 



sec 




dX 1 


( 11 ) 


For ozone, the effective absorptance a(O^) over the entire solar spectrum 
is given by Lacis and Hansen (1974) as 


a(03> = 


0.02118U 


1.082U 


1 + 0.042U + 3.23 X lo"'^ u^ 


(1 + 138. 6u) 


0.805 


0. 0658u 


1 + (103. 6u) 


(12) 


where 


u = £ sec 0 (13) 

The effective absorptance for water vapor according to Lacis and Hansen 
(1974) is 




2.9w 

(1 + 141.5w)°'^^^ + 5.925W 


(14) 


where 


w = X cos 0 
w 


(15) 
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Ozone absorption and Rayleigh scattering occur in the same region of the 
spectrum, but ozone absorption is at high altitudes where there is little 
Rayleigh scattering. Therefore the solar radiation reaching the surface if 
there were only attenuation by ozone and Rayleigh scattering would be 


F(R,H 20 ) = Iq cos 0 exp 




T0(R) 


1013 


sec 


e [l - a( 0 ^)] 


(16) 


Water vapor absorption occurs in the near-infrared region of the spectriam where 
ozone absorption and Rayleigh scattering are negligible. Thus the attenuation 
by water vapor can be subtracted directly from equation (16) to yield the direct 
solar beam of radiation at the surface as 


I 



P 

1013 


sec 


0 


[l - a( 03 )] 


a(H20) 


(17) 


The diffuse radiation arriving at the surface consists of the solar radia- 
tion scattered downward by the molecules and aerosols. With scattering by aero- 
sols neglected, scattering downward to the surface can be approximated by the 
following equation (Paltridge and Platt, 1976) ; 

0 

D = I^Iq - I) COS 0 - Iq cos 0 aj^(0)J 

+ I cos 0 a (18) 

g R 

In equation (18), (Iq - I) is the radiation removed from the direct beam by 
scattering, Iq cos 0 aj^(0)- is the upward scattered radiation, and Uj^(0) is 
an approximate coefficient for upward scattering (i.e. atmospheric reflection) 
given by Lacis and Hansen (1974) as 


(l + ag4) 


aR(0) 


,0.28 

1 + 6.43 cos 0 


(19) 


— * 

The term Oi„o,„ takes into account the downward scattering of the diffuse radia- 
tion reflected from the surface; is the ground albedo, and the 

reflection coefficient for downward scattering, about 0.0685. The last term in 
equation (18) is the direct beam radiation which is reflected from the Earth 
and then scattered back downward. 
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EMPIRICAL MODELS 


Clear Sky: Zenith Angle Variation 

Under clear sky conditions, zenith angle is the main factor affecting the 
solar radiation arriving at the surface at a given location. In this case 
empirical relations can represent the long time average of direct, diffuse, and 
global solar radiation with good accuracy. For example, Paltridge and Platt 
(1976) derived empirical functions from data taken at Aspendale, Victoria, 
Australia, for the 5-year period from 1967 to 1972. The data and empirical 
curves are shown in figures 3 to 5. The empirical functions are as follows. 

For global solar radiation on a horizontal surface, at Aspendale 
(latitude 38°02's), the hourly total (mW-hr/cm^) is 


F = 1.0 + 141.1 sin Y - 31.0 (sin Y) 


1/2 


( 20 ) 


where Y is the solar elevation angle. The direct normal solar radiation is 
given by 


I = 100. 0[l - exp(-0.06Y)] (21) 

and the diffuse radiation on a horizontal plane by 

D = 0.5 + 9.6[l - exp(-0.05Y)] (22) 


Cloudy Skies 

It is difficult to treat clouds in solar energy calculations. Detailed 
studies require knowledge of cloud drop size distribution and cloud total water 
content. Although clouds generally diminish the solar radiation incident on the 
surface, the surface radiation may actually be increased under partially cloudy 
conditions if the solar disk is not occulted. 

Table V indicates the percent decrease in fliax of total surface radiation 
due to continuous cloud cover of various types as a function of the solar 
zenith angle (Kondrat'ev, 1973). 

An empirical expression that yields reasonably accurate results for mean 
annual and mean monthly studies is (Berliand, 1961) for global radiation: 

F = Fq( 1 - an - bn^) (23) 
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where is the clear sky radiation, b = 0.38, a is a coefficient that 

varies with latitude as shown in table VI, and n is cloud cover in tenths. 
This relation should be used only for long time averages and not for calcula- 
tions on a daily basis. 


Models of the Air Resources Laboratory (NOAA) 

The Air Resources Laboratory of NOAA has developed empirical models of 
global solar radiation on a horizontal plane. The models were developed for the 
rehabilitation of historical data archives for 26 stations in the NOAA solar 
radiation network (Solmet, 1978) . The data tapes consisted of hourly records 
of observed weather data and 

1. Integrated hourly global solar radiation 

2. Duration of sunshine in minutes 

3. Cloud opacity in tenths 

4. Cloud cover in tenths 

5. Sky condition observations 

Five regression equations were developed and used for those stations with an 
appropriate data set for the regression analysis. The form of the equations 
for each model and examples of the parameters are as follows: 

1. Clear Sky 

o 3 

SRC = Aq + A^ cos(ZA) + A^ cos (ZA) + A^ COS (ZA) (24) 

2 

where SRC is the clear sky hourly integral global solar radiation (kJ/m ) 
and ZA is zenith angle. The coefficient Aq is determined for mornings 
and afternoons for each month of the year; and Aj^, A 2 , and A^ are esti- 
mated separately for mornings and afternoons, one set of each for the year. 

2. Sunshine and Opaque Cloud (preferred estimating function) 

SR = SRC^Bq + Bj^(SS) + B2(0PQ) + B 3 (OPQ) ^ + B^ (OPQ) ^ + B^ (RN^ (25) 

where SR is the hourly global solar radiation (kJ/m ) , SS is the number 
of minutes of sunshine divided by 60, OPQ is cloud opacity (0.1, 0.2, etc.), 
and RN is a precipitation variable which is 0 for no precipitation and 1 
when precipitation is reported. 
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3 . Opaque Cloud Only (Scime as 2 without sunshine term) 


SR = SRc|cq + C 2 ( 0 PQ) + C 3 ( 0 PQ)^ + (OPQ) ^ (RN)j 

4 . Sunshine Only (converse of 3 ) 


( 26 ) 


SR = SRC[Dq + D^(SS) + D^CRN)] 


5 . Sky Condition 


SR = SRC 


Eq + ^ Ei(SC^) + Eg(RN) 

j=l 


( 27 ) 


( 28 ) 


where each SC^^ is 0 or 1 depending on whether the sky cover listed below is 
present at any one of up to four levels: 

SCj^ thin scattered; 0.1-0 . 5 cover 

SC2 opaque scattered; 0.1-0 . 5 cover 

SC3 thin broken; 0.6-0 . 9 cover 

SC^ opaque broken; 0.6-0 . 9 cover 

SC3 thin overcast; 1.0 cover 

SCg opaque overcast; 1.0 cover 

SC-^ partial or total obscuration 

The coefficients for the regression equations are available on punched cards; 
an example is listed below. The punched cards are arranged to provide a maximum 
of 12 fields of length 6 (including the decimal point) . The last 8 columns of 
each card are identification. There are 7 cards for each station as follows: 
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Card 

order 

Parameter 

Niimber of 
coefficients 

Number of punched 
fields on card 

Description 

1 

^0 


12 

12 

Clear sky; Jan-Dec, 





morning 

2 

^0 


12 

12 

Clear sky; Jan-Dec, 
afternoon 

3 

^l'^2' 

^3 

6 

6 

Clear sky; first 3 are 
morning; last 3 are 
afternoon 

*4 

Co - 

C5 

5 

6 

Opaque cloud only: 

SS coef. = 0.0 (CjL = 0) 

*5 

^0 - 

^5 

6 

6 

Sunshine and opaque 
cloud coef. 

*6 

- 


3 

6 

Sunshine only: opaque 

cloud coefs. = 0.0 
(d2 = d^ = d^ = 0) 

7 

^0 - 

^8 

1 

9 

9 

Sky cond. coef. 


*Cards 4, 5, and 6 have common ordered positions for the constant term, 
SS, OPQ, OPQ^, OPQ^, and RN. For example, on card 4 the sunshine coefficient 
(field 2) is always zero. 


An example of a set of coefficients is 
Washington, D.C. 



12. 

1. -30. 

-89. 

-180. -219. -192. -128. 

-65. 

-37. 

-26. 

WSH 

CLRM 

9. 

20. -30. 

-87. 

-172. -176. -195. -131. 

-67. 

-59. 

-15. 

WSH 

CLRA 

2810. 

-1540. 2100. 

3570. 

-2060. 




WSH 

CLRZ 

0.000 

-.474 0.920 

1.070 

-.214 




WSH 

0 

0.000 

0.000 0.000 

0.000 

0.000 




WSH 

SS.O 

0.000 

0.000 0.000 

0.000 

0.000 




WSH 

SS 

0.048 

-.005 -.007 

-.141 

-.055 -.471 -.131 -.251 




WSH 

SKYC 








THE INFERENCE OF SOLAR RADIATION AT THE SURFACE FROM SATELLITE DATA 


The technique involved in the application of satellite images to the pre- 
diction of solar radiation arriving at the Earth's surface has three steps: 

1. Establish a set of equations which can be used to relate the data from 

the satellite images to the insolation. It may also be necessary to 
include some conventional meteorological parameters in the equations . 

2. Use satellite data for a period of time plus surface solar energy 

measurements taken at the same time to determine unknown constants in 
these equations . 

3. Then use satellite image data with the equations to predict the solar 

radiation incident upon the surface. 

In the following, one of the more recent attempts, estimating hourly and daily 
global solar radiation from GOES geostationary satellite images, is discussed. 


The Great Plains Experiment 

The National Environmental Satellite Service and the Great Plains Agri- 
cultural Council have carried out an experiment to determine incoming solar 
radiation in the U.S. Great Plains area from geostationary satellite images 
(Tarpley, 1979) . 

The region selected for the experiment, shown in figure 6, was bounded by 
two latitude circles, 29° and 49° N, and two meridians of longitude, 95° and 
105° W. Twenty-two sites were instrumented with pyranometers for measurement 
of incident global surface energy at the points shown in the figure. 

The satellite data used were visible (0.55 to 0.75 ym) images from the 
GOES visible and infrared spin scan radiometer (VISSR) , having 8 km resolution 
at Nadir. Relative brightness levels were digitized (0 to 63 counts). Hourly 
images were used for 7 to 10 hours per day distributed about local noon for the 
period from June 7, 1977, to August 15, 1977. 

Surface meteorological data used were surface pressures and total pre- 
cipitable water. 

The grid used for analyses of the satellite images is shown in figure 7. 
The Great Plains area was divided into targets approximately 50 km on a side, 
each made up of 7 x 6 arrays of 8 km pixels. 

Quantities used as variables for determination of surface insolation were 

1. Surface pressure, a measure of total air mass, which therefore should 
be correlated with the depletion of solar radiation by molecular 
scattering 
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2. Precipi table water which provides information on absorption and 

scattering by water vapor 

3. Mean target brightness which is a measure of reflected and scattered 

radiation due to all causes 

4. Cloud amount and brightness which indicate the amount of radiation 

reflected by the clouds. 


Determination of Cloud Free Brightness Levels 

An initial task was to determine cloud free brightness in the images, by 
regression against functions of solar zenith angle X and the azimuth angle (p 
between the Sun and the satellite. Cloud free brightness was calculated on a 
2° grid with interpolation to determine values at 0.5° target intervals. For 
this purpose, 100 observations per target were used for 27 days in May 1977, 
just previous to the period of the test. 

In order to obtain the regression coefficients, cloud contaminated observa- 
tions were eliminated by an automatic cloud detection and elimination procedure, 
along with fitting the data to the cloud free regression equation of the form 

2 

B = A + b cos X + c sin X cos sin X cos <() (29) 

The cyclic procedure, shown in figure 8, produced cloud free data sets and 
regression coefficients after two iterations of the procedure, with the final 
standard error of the estimate for the regression equation of 0.5 to 0.9 counts. 


Determination of Cloud Parameters and Target Brightness 

The cloudiness of each pixel was established according to the following 
criteria: it was clear if the brightness was less than or equal to B(X,4>) + 3, 

partly cloudy if the brightness was greater than B(X,({)) + 3 but less than or 
equal to B(X,<t>) + 5, and cloudy if the brightness was greater than B(X,<J>) + 5. 

The cloud fraction for the target was then calculated from 


O.5N2 + 

^1 + t N3 


(30) 


where N 2 , and were the number of pixels in the clear, partly cloudy, 

and cloudy cases. 

The mean cloud brightness I^ld taken to be the mean of the brightness 

for the cloudy cases. The quantity 


143 


(31) 


Bq = B (45°, 105°) 


was also calculated as the normalized clear brightness. 


Estimation of Hourly Insolation 

The form of the regression equations for hourly insolation was suggested by 
consideration of the equation of conservation of radiant energy 

~ 2a 2g 

where the subscripts indicate the following: 

0 at the top of the atmosphere 

R reflected back into space 

A absorbed in the atmosphere 

G absorbed at the ground 


Qg = Qgd - a) 

where Qg is the insolation at the surface and ex is the surface albedo. 
Substituting from equation (32) for Qq results in 


“ 1 - 2r - Qa> 


The regression equations used were determined by cloud amount; they were 




Qc = a + b cos X cV + dn + et-:^ 


(n < 0.4) 


Qo = a + b cos X + cn 


(0.4 ^ n < 1.0) 


Qg = a + b cos X + c 
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where 


2s 

X 

n 


'•cld 


'0 


hourly surface insolation, ly (1 langley = 41.84 kJ/m^) 

local solar zenith angle 

transmittance of the clear atmosphere 

fractional cloud amount 

mean target brightness 

mean cloud brightness 

predicted clear brightness 

normalized clear brightness, B (45°, 105°) 


a,b,c,d,e the regression coefficients 

Precipitable water and surface altitude Z were used in the determination 
of the clear target transmission: 


= lb \l) Ip 
^ ^ws^wa^r 


(36) 


where the transmission due to water vapor scattering, water vapor absorption, 
and Rayleigh scattering, respectively, were taken to be; 


^ws 


1 - 0.00225wm 


= 1 - 0.077 (wm) 


0.3 


ip^ = 0.972 - 0.0862m + 0.00933m^ 


(37) 

(38) 

(39) 


where w is the precipitable water and m, the optical air mass, is given by 


m = sec X exp 


(- 




8243/ 


(40) 


Results 

To determine the technique for estimating surface insolation from cloud 
images, all of the data were included. This is referred to as the dependent 
development data set. In order to estimate how accurately insolation could be 
estimated under operational conditions, the data were divided into two parts, 
the dependent set with only 4 pyranometer sites and the independent set with 
18 pyranometers to verify the regression equation. 
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The values of the regression coefficients and their standard errors are 
listed in table VII for the dependent developmental data set- For this set, 
satellite estimates from the regression equations are plotted versus actual 
measurements in figures 9(a), 9(b), and 9(c) for n < 0.4, 0.4 - n < 1.0, and 

n = 1.0, respectively. Correlation coefficients and standard errors of the 
estimates are shown in the figures. 

Daily cumulative insolation was also estimated, using the zenith angle to 
interpolate and fill in gaps in the data. These results are shown in figure 10. 

Finally the statistics of estimated total insolation for four different 
cases are shown in table VIII. Cases considered are 

1. Seven or more pictures per day, dependent data 

2. Seven or more pictures per day, independent data 

3. Two pictures per day, dependent data 

4- One picture per day, dependent data 

In each of the first three cases, the standard errors are less than 10 percent 
of the mean. In the last case, with only one picture per day, the standard 
error was about 20 percent of the mean. In this last case, sampling was 
insufficient. 

The most significant problem with the technique is the overestimation of 
surface insolation under cloudy conditions (see fig. 9(c)). 
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PAGE iii 

OF POOR QUALITY 


TABLE I.- SPECTRAL ABSORPTION COEFFICIENTS OF OZONE (cm 

[From Leckner, 1978*] 


X 

k 

X 

k 

X 

k 

X 

k 

X 

k 

0.290 

38.000 

0.445 

0.003 

0.515 

0.045 

0.585 

0.118 

0.700 

0.023 

0.295 

20.000 

0.450 

0.003 

0.520 

0.048 

0.590 

0.115 

0.710 

0.018 

0.300 

10.000 

0.455 

0.004 

0.525 

0.057 

0.595 

0.120 

0.720 

0.014 

0.305 

4.800 

0.460 

0.006 

0.530 

0.063 

0.600 

0.125 

0.730 

0.011 

0.310 

2.700 

0.465 

0.008 

0.535 

0.070 

0.605 

0.130 

0.740 

0.010 

0.315 

1.350 

0.470 

0.009 

0.540 

0.075 

0.610 

0.120 

0.750 

0.009 

0.320 

0.800 

0.475 

0.012 

0.545 

0.080 

0.620 

0.105 

0.760 

0.007 

0.325 

0.380 

0.480 

0.014 

0.550 

0.085 

0.630 

0.090 

0.770 

0.004 

0.330 

0.160 

0.485 

0.017 

0.555 

0.095 

0.64 0 

0.079 

0.780 

0.0 

0.335 

0.075 

0.490 

0.021 

0.560 

0.103 

0.650 

0,067 

0.790 

0.0 

0.340 

0.040 

0.495 

0.025 

0.565 

0.110 

0.660 

0.057 

0.800 

0.0 

0.345 

0.019 

0.500 

0.030 

0.570 

0,120 

0.670 

0.048 

0.810 

0.0 

0.350 

0.007 

0.505 

0.035 

0.575 

0.122 

0.680 

0.036 

0.820 

0.0 

0.355 

0.0 

0.510 

0.040 

0.580 

0.120 

0.690 

0.028 

0.830 

0.0 


TABLE II.- EFFECTIVE ABSORPTION COEFFICIENTS OF WATER VAPOR (cm^g 


[From Leckner, 1978*] 


X 

k 

w 

X 

k 

w 

X 

k 

w 

X 

k 

w 

X 

k 

w 

0.69 

0 . 160F-01 

0.64 

0 . 1558+00 

0.99 

0 . 1258+00 

1.70 

0 . 5108+00 

2.90 

0 . 6508+03 

0.70 

0 . 2401-01 

0.85 

0 . 3008-02 

1.00 

0 . 2508-02 

1,75 

0 . 4008+01 

3.00 

0 . 2408+03 

0.71 

6 . 1258-01 

0.86 

0 . 1008-04 

1.05 

0 . 1008-04 

1 .80 

0 . 1308+03 

3.10 

0 . 2308+03 

0.72 

0 . 1008+01 

0.87 

0 . 1008-04 

1.10 

0 . 3208+01 

1.65 

0 . 2208+04 

3.20 

0 . 1008+03 

0.73 

0 . 8708+00 

0.88 

0 . 2608-02 

1.15 

0 . 2308+02 

1 .90 

0 . 1408+04 

3.30 

0 . 1208+03 

0,74 

0 . 6108-01 

0.89 

0 . 6308-01 

1 .20 

0 . 1608-01 

1.95 

0 . 1608+03 

3.40 

0 . 1958 + 02 

0.75 

0 . 100F-02 

0.90 

0 . 2108+01 

1.25 

0 . 1808-03 

2.00 

0 . 2908+01 

3.50 

0 . 3608+01 

0.76 

0 . 1008-04 

0.91 

0 . 1608+01 

i.30 

0 . 2908+01 

2.10 

0 . 2208+00 

3.60 

0 . 3108+01 

0.77 

0 . 1008-04 

0.92 

0 . 1258+01 

1.35 

0 . 2008+03 

2.20 

0 . 3308+00 

3.70 

0 . 2508+01 

0.78 

0 . 6008-03 

0.93 

0 . 2708+02 

1 .40 

0 . 1 108+04 

2.30 

0 . 5908+00 

3.60 

0 . 1408+01 

0.79 

0 . 1758-01 

0.94 

0 . 3008+02 

1 .45 

0 . 1508+03 

2.40 

0 . 2038+02 

3.90 

0 . 1708+00 

0.80 

0 . 3608-01 

0.95 

0 . 4108+02 

1.50 

0 . 1508+02 

2.50 

0 . 3108+03 

4.00 

0 . 4508-02 

0.81 

0 . 3308+00 

0.96 

0 . 2608+02 

1.55 

0 . 1708-02 

2.60 

0 . 1508+05 



0.82 

0 . 1538+01 

0,97 

0 . 3108+01 

1.60 

0 . 1008-04 

2.70 

0 . 2208+05 



0.83 

0 . 6608+00 

0.98 

0 . 1468+01 

1 .65 

0 . 1008-01 

2.60 

• 0 . 6008+04 




TABLE III.- TYPICAL ANGSTROM TURBIDITY COEFFICIENTS (6q) 

The original data are converted from Schuepp's turbidity coefficient B 
through the relation 3 q = 0.935B. The three decimal digits in the 
table are a result of this conversion and do not indicate the 
accuracy; note that 6 = from Leckner, 1978* 


Latitude 

60°N 

45°N 

30°N 

0° 

Low limit 

0.01 

0.047 

0.047 

0.047 

Medium 

0.057 

0.077 

0, 097 

0.117 

High limit 

0.093 

0.187 

0.375 

0.375 


*Reprinted with permission from Sol. 
Pergamon Press, Ltd. 
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TABLE V.- RELATIVE DECREASE (IN PERCENT) OF THE FLUX OF TOTAL RADIATION 
FOR THE CASE OF A CONTINUOUS CLOUD COVER OF DIFFERENT TYPES AND 
DIFFERENT SOLAR ALTITUDES IN COMPARISON WITH THE 
CORRESPONDING VALUES FOR A CLOUDLESS SKY 
[ From Kondrat ' ev , 1973 ] 


Cloud 

type 

Solar altitude, deg 

5 

10 

20 

30 

40 

50 

Ci 

44 

50 

34 

22 

10 

2 

Cs 

33 

39 

51 

35 

20 

10 

Ac 

33 

39 

46 

55 

46 

38 

As 

44 

50 

59 

64 

63 

63 

St fr . 

78 

83 

80 

77 

77 

76 

Sc 

89 

78 

68 

71 

72 

73 

St 

78 

78 

80 

81 

83 

84 


TABLE VI.- COEFFICIENT a IN EMPIRICAL EQUATION FOR MEAN MONTHLY AND MEAN 
ANNUAL GLOBAL RADIATION UNDER CLOUDY SKIES 


Lat. 

0° 

5° 

o 

o 

15° 

20° 

25° 

0 

o 

n 

35° 

0 

o 

a 

0.38 

0.40 

0.40 

0.39 

0.37 

0.35 

0.36 

0.38 

0.38 

Lat. 

45° 

50° 

55° 

60° 

65° 

70° 

75° 

80° 

85° 

a 

0.38 

0.40 

0.41 

0.36 

0.25 

0.18 

0.16 

0.15 

0.14 
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TABLE VII.- REGRESSION COEFFICIENTS, THEIR STANDARD ERRORS, AND 
OTHER STATISTICS FOR HOURLY INSOLATION REGRESSIONS 
[From Tarpley, 1979] 



a, ly 

b, ly 

Sb- ly 

c, ly 

Sc» ly 

d, ly 

Sd' ly 

e, ly 

Sg/ ly 

Multiple- correlation 
coefficient . . . 


Standard error of 
the estimate . . 


Number of cases 


Clear cases 
(n < 0.4) 

Partly cloudy 
(0.4 ^ n < 1.0) 

Cloudy 
(n = 1.0) 

-19.33 

-9.57 

-6.56 

87.08 

94.54 

87.68 

0.71 

1.78 

3.16 

27.58 

-7.62 

-7.50 

2.33 

0.28 

0.35 

-10.48 



1.33 



-6.37 



0.79 



0.94 

0.77 

0.70 

5.56 

12.19 

11.36 

5736 

2127 

822 


TABLE VIII.- STATISTICS OF ESTIMATED TOTAL DAILY INSOLATION 

[From Tarpley, 1979] 



Mean error, ly .... 
Correlation coefficient 
Standard error, ly . . 
Mean daily insolation 
(pyranometer) , ly . . 
Number of cases .... 


Seven 
or more 
pictures 
per day 
(dependent 
data) 

Seven 
or more 
pictures 
per day 
(independent 
data) 

Two 

pictures 
per day 
(dependent 
data) 

10.7 

17.5 

12.4 

0.91 

0.90 

0.88 

50.4 

53.5 

58.5 

586.2 

594.4 

586.7 

896 

721 

1036 


One 

picture 
per day 
(depender 
data) 



5 





















OF POOR QUALITY 



Figure 1.- Calculated direct irradiance on a horizontal surface at air mass 2 
and two values of turbidity coefficient B, compared with measured values. 
The missing weak bands in the calculated curves are indicated with arrows. 
(From Leckner, 1978; reprinted with permission from Sol. Energy, vol. 20, 
copyright 1978, Pergamon Press, Ltd.) 



Wavelength um 


Figure 2.- Calculated global and diffuse irradiances on a horizontal surface 
at air mass 2 and two values of turbidity coefficient B» compared with 
measured values. (From Leckner, 1978; reprinted with permission from Sol. 
Energy, vol. 20, copyright 1978, Pergamon Press, Ltd.) 
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Figure 3.- Global solar radiation (hourly total) in 
clear skies as a function of solar elevation Y 
at Aspendale for the period from 1967 to 1972. 
(From Paltridge and Platt, 1976.) 
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Figure 4.- Direct radiation (hourly total) in clear skies as a function of 
solar elevation y Aspendale for the period from 1967 to 1972. 

(From Paltridge and Platt, 1976.) 


40 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



Solar elevation (deg) 

Figure 5.- Diffuse radiation (hourly total) in clear skies as a function of 
solar elevation angle y at Aspendale for the period from 1967 to 1972. 
Note the definition of "clear sky" in the text. (From Paltridge and 
Platt, 1976.) 
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ORIGINAL PAGE iS 
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Figure 6.- The part of the Great Plains where satellite data 
were collected (dark line) . The pyranometer sites used 
in this study are marked with dark circles. (From 
Tarpley, 1979.) 



Figure 7.- Grid used for analysis of satellite images. 
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Lgure 8.- Flow diagram for processing target brightness 
data into regression coefficients. (From Tarpley, 
1979 . ) 
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(a) 5736 clear cases (b) 2217 partly cloudy cases 

(n < 0.4) , (0.4 ^ n < 1.0) . 



(c) 822 totally cloudy cases (n = 1.0). 


Figure 9.- Hourly insolation measured by the 22 pyranometers plotted against 
satellite estimates (dependent data) . If there were no errors, all points 
would lie on the diagonal line. Correlation coefficients r and standard 
errors 0 are shown. (From Tarpley, 1979.) 
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Figure 10.- Daily insolation measured by the 22 pyranometers 
(dependent data) plotted against satellite estimates 
(^7 pictures per day) . Diagonal line shows perfect agree 
ment between pyranometers and satellite estimates. 
(Correlation coefficient r and standard error a are 
shown. (From Tarpley, 1979.) 




